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ABSTRACT 


The  velocity  structure  in  the  earth's  mantle  is  de¬ 
termined  using  the  trave±  times  and  travel -time  slopes  of 
P-waves  recorded  at  LASA.  The  results  show  anomalous 
velocity  gradients  at  depths  of  700,  1200,  and  1900  km. 
Both  surface  wave  dispersion  data  and  dt/dA  data  show  the 
presence  of  lateral  heterogeneities  in  the  upper  and 
lower  mantle. 

Theoretical  studies  of  acoustic  and  gravity  wave 
propagation  in  the  atmosphere -ocean  systems,  tfffect  of 
aftershocks  on  free  oscillations,  and  response  of  layered 
spherical  ^arth  to  point  sources  have  been  described. 
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I.  INTRODUCTION 

This  is  the  first  Annual  Report  summarizing  cur 
work  in  "Research  in  Seismology"  under  AFOSR  contract. 

During  the  year  we  have  initiated  the  work  on 
a  number  of  problems.  A  few  of  these  have  been  completed, 
but  most  of  them  are  still  in  progress.  In  this  report  we 
describe  briefly  those  projects  which  have  been  advanced 
to  the  stage  where,  at  least,  some  preliminary  results 
are  available. 

Our  efforts  have  been  directed  to  several  areas; 
These  are: 1)  utilization  of  data  from  LaSA  and  other 
sources  to  study  the  structure  of  the  earth's  interior, 
l)  theoretical  studies  of  seismic  wave  propagation  and 
detection  methods,  3)  source  mechanism  studies,  A)  estab¬ 
lishment  of  a  reference  seismic  station  equipped  with  long 
period  instruments.  In  the  following  chapters  each  one 
of  these  areas  will  be  described  briefly. 
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II.  STRUCTURE  OF  THE  EARTH’S  MANTLE 
A .  Body  wave  results. 

Seismic  body  and  surface  waves  generated  by 
earthquakes  and  underground  explosions  have  been  utilized 
to  study  the  structure  of  the  earth's  mantle.  Phase  and 
group  velocities  of  the  Rayleigh  and  Love  waves  in  the 
period  range  of  80-600  seconds  were  utilized  to  interpret 
the  shear  wave  structure  in  the  upper  mantle.  The  travel 
times  of  P-waves  from  explosions  and  the  travel  time 

I  'T 

slopes  )  measured  from  LASA  recordings  were  util¬ 

ized  to  obtain  the  compressional  velocity  profile  in  the 

/ 

mantle . 

Because  of  its  large  aperture,  geometry,  and 
matched  recording  system^.  LASA  provides  an  excellent  means 
of  measuring  the  travel  time  slopes  (  “VdA  )  directly. 

This  in  turn  can  be  interpreted  in  terms  of  velocity  dis¬ 
tribution.  We  have  carried  out  such  a  study  utilizing 
the  data  computed  from  LASA  recordings  of  events  at 

p 

distances  greater  than  A  ~  "cO  ,  together  with  the  travel 
times  from  the  Long  Shot  explosion.  To  minimize  the  effect 
of  lateral  inhomogeneities  in  the  mantle,  the  analysis 
was  confined  to  events  arriving  from  a  narrow  wedge  of 
azimuth  300-320°,  measured  clockwise  from  the  north  at 
LASA.  Figure  1.  defines  the  bard  which  this  wedge  covers. 
Because  of  the  extreme  seismicity  of  the  region^  a  majority 
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of  teleseismir  events  detected  at  LASA  fall  in  this  range. 
Furthermore,  the  Long  Shot  explosion  provides  accurate 
travel  time  data  for  the  interpretation,  and  for  making 
static  correction  to  curve,  inis  ambiguity  arises 

because  of  the  station  corrections  at  LASA. 

The  measured  are  compared  with  the  Jef- 

freys-Bullen  model  in  Figure  2.  The  discrepencies  are 
well  defined  and  obvious.  A  new  velocity  model  was  com¬ 
puted  for  the  mantle  to  fit  the  measured  and  the 

travel  time  residuals  from  the  Long  Shot.  The  observed 
residuals  and  the  curve  computed  using  the  new  velocity 
model  for  the  mantle  is  given  in  Figure  4.  The  portion 
of  the  model  between  the  surface  and  the  depth  of  600  km 
is  based  on  the  surface  wave  data  and  other  body  wave 
studies.  The  portion  below  that,  however,  satisfy  both 
the  travel  time  residuals  and  measurements.  Since 

the  latter  is  especially  sensitive  to  velocity  gradients 
at  depth,  the  changes  in  velocity  gradients  shown  at  depths 
of  700,  1200,  and  1900  km  are  reliable  and  the  most  inter¬ 
esting  features  of  the  lower  mantle  structure.  These  and 
the  anomalous  region  at  350  km  depth,  constitute  a  set  of 
"discontinuities"  for  P-wave  velocity  profiles  in  the  mantle. 

The  indications  are  such  that  there  are  lateral 
variations  in  the  velocity  structure  of  the  lower  mantle. 

In  Figure  5  /&t>  vs*  A  plots  for  arrivals  to  LASA 
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from  two  different  azimuths  are  shown.  The  0  "300-320° 
azimuth  corresponds  to  the  Aleutian-Kuriie  paths  outlined 
in  Figure  and  the  0  “120-140°  azimuth  to  Mexico - 
Central  and  South  America.  The  differences  in  the  slopes 
of  vs*  ^  plots  are  especially  outstanding  at 

regions  A  ^50°  >  indicating  lateral  variations  deep  in 
the  lower  mantle. 


B.  Lateral  Inhomogeneities  in  the  Upper  Mantle  from 
Surface  Waves. 

The  lateral  inhomogeneities  in  the  uppermost 
portion  of  the  mantle  have  been  established  through  seis¬ 
mic  refraction  studies.  At  greater  depths,  complications 
arise  in  the  use  of  travel  times  because  of  low  velocity 
zone.  There  are  some  regional  disoersion  data  and  mantle 
structures.  However,  most  data  at  periods  longer  than 
100  seconds  are  computed  ovei  mixed  paths  including  con¬ 
tinental  shields,  oceans,  and  tect-onic  areas.  With  certain 
assumptions,  pure  path  phase  velocities  have  been  extracted 
from  composite  Love  wave  data  and  the  corresponding  struc¬ 
tures  computed.  These  are  shown  in  Figures  6  and  7.  The 
differences  in  the  upper  mantle  shear  velocity  profiles 
are  outstanding.  Oceanic  areas  seem  to  have  the  lowest 
velocities  in  the  low  velocity  zone  and  the  continental 
shields  the  highest  tectonic  regions  fall  in  between. 


A  technique  has  been  utilized  to  estimate  the 
temperature  differences  between  oceans  and  continental 
shields  from  the  velocity  differences  shewn  abeve.  With 
simplifying  assumptions  this  can  be  written  as, 


(^-1 

A  /  /  /  -j— 

where  Av  ■-  V£-V0  is  velocity  difference,  A  I  (2  ^tem¬ 
perature  coefficient  of  velocity.  The  last  quantity  must 
be  determined  experimentally.  The  temperature  differences 
were  computed  taking  various  values  for  (  )p within 

the  limits  of  available  laboratory  data.  These  are  shown 
in  Figure  8„  The  oceanic  upper  mantle  is  characterized 
by  higher  temperatures,  with  the  maximum  temperature  dif¬ 
ference  being  at  about  110  km  depth. 


W  Affc-rv+j  d£) 
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III.  THEORETICAL  STUDIES. 

A.  Pulse  Propagation  in  Coupled  Systems 

The  theoretical  and  numerical  techniques  of 
pulse  propagation  in  a  laterally  homogeneous  layered  at¬ 
mosphere  is  extended  to  coupled  systems  such  as  the  one 
formed  by  the  atmosphere  and  the  ocean.  The  theory  was 
applied  to  the  air-sea  waves  excited  by  the  explosion 
of  the  volcano  Kiakatoa. 

In  numerical  calculations  a  standard  ARDS  at¬ 
mosphere  underlain,  by  a  constant-density,  constant-velocity 
ocean  was  assumed.  The  atmosphere  was  terminated  by  a 
free  surface  and  the  ocean  bottom  was  assumed  to  be  per¬ 
fectly  rigid.  The  explosive  source  was  placed  at  the  sur¬ 
face  of  the  ocean. 

The  dispersion  curves  for  the  gravity  (GR)  and 
acoustic  (S)  waves  in  the  atmosphere  and  the  classical 
gravity  surface  waves  in  the  ocean  (GW0)  were  computed. 
These  are  shown  in  Figure  9  and  10.  The  results  indicate 
that  free  waves  in  the  atmosphere  do  exist  with  phase 
velocities  near  those  of  GWQ  (C  at  long  periods). 

This  leads  to  an  efficient  transfer  of  energy  from  the 
atmosphere  to  the  ocean  and  produces  sea  waves  with  am¬ 
plitudes  several  times  larger  than  hydrostatic.  The 
phase  velocity  curves  for  an  atmosphere -ocean  coupled 
system  are  shown  in  Figure  11.  The  modif ications  from 
that  of  an  atmosphere  alone  can  be  seen  from  the  compari¬ 
son  of  Figures  9  and  11. 
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Synthetic  barogram  and  a  synthetic  marigram  were 
constructed  using  a  surface  source  with  a  time  function 
of  a  single-cycle  sine  wave  of  period  40  minutes.  These 
records  ere  shown  in  Figure  12.  The  distance  corresponds 
from  Krakatoa  to  San  Francisco.  The  initial  pulse  is  a 
transient  visible  both  as  an  atmospheric  and  oceanic  dis¬ 
turbance  of  sea  level  with  arrival  of  GW0  mode  with  hardly 
any  pressure  pulse  accompanying  it. 

The  scaling  laws  of  atmospheric  nuclear  explo¬ 
sions  applied  to  the  Krakat-oa  shows  that  the  latter  was 
equivalent  to  10U-150  megaton  explosion. 


References 


Harkrider,  D.  G.  and  F.  Press.  The  Krakatoa  air-sea 
waves:  An  example  of  pulse  propagation  in  coupled 
systems,  Geophysical  Jour..  1966  (in  press). 

Press,  F.  and  D.  G.  Harkrider.  The  air-sea  waves  from 
the  explosion  of  Krakatoa,  Science.  1966,  (in  press). 


B.  Effect  of  Aftershocks  on  the  Measured  Period  and  Q  of 
Free  Oscillations. 

In  the  estimations  of  the  Q  and  the  eigen- 
frequencies  of  the  free  oscillations  of  the  earth  using 
power  spectra,  the  effects  of  aftershocks  were  not  taken 
into  account.  Since  most  great  earthquakes  are  followed 


by  a  large  number  of  aftershocks  of  varying  magnitudes, 
it  would  be  desirable  to  estimate  the  effects  of  these 
aftershocks  on  the  measured  free  oscillation  spectra. 

In  formulating  the  problem  of  the  main  event 
and  the  aftershocks  as  a  series  of  K  events  each  of 
which  occurs  at  time4 

-=Za(-V>  f  U-V) 

Let  g(w)  e  be  the  spectrum  of  f  (  )  and^be 

t“h 

frequency  independent  amplitude  factor  for  K  event. 

Now,  let's  compute  the  spectral  density  of  the  time  series 

1  G(aO(  }  in  terms  of  I  '  l  .  Taking  the  Fourier 

transform  and  squaring  we  obtain, 

/->  2.  <2  -v  ~  LU* 

Gc^Go^r=  cp;|  2.  e  ^ 

If  it  is  assumed  that  both  the  amplitude  and  probability 
of  the  occurence  of  events  with  decreases  exponent- 


ially  as 

,  .  A  -<Ta 

(X and 

,  the  result  reduces 

to, 

[  GTo; 

r= 

1KaMk"  + 

1C  (K  -  i)  €'*’ 

*0  ' 

(cr-^%  u)* 

This  shows  that  |£i(.td)|  differs  from  [^Cwj| 

Taking  the  actual  aftershock  distribution  from 
Chile  Earthquake  of  1960,  and  Alaska  Earthquake  of  1964, 
the  effects  of  aftershocks  on  the  Q  and  eigenfrequency  of 
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main  shock  were  computed.  The  results  of  one  mode  are 
shown  in  Figure  13,  where  a  shift  in  peak  frequency  of 
about  0.01  minutes  can  be  seen.  The  Q  is  increased  i  om 
205  to  250. 

On  the  basis  of  calculations  it  is  concluded 
tha  •  the  aftershock  sequences  can  affect  the  position 
and  shape  of  spectr'l  peaks  by  a  finite  amount.  The  mag¬ 
nitude  of  the  change  of  the  peak  frequency  is,  in  general, 
too  small  to  affect  the  interpretation.  The  errors  in¬ 
troduced  in  Q  measurements  are  greater  and  become  sig¬ 
nificant  at  shorter  periods. 
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C .  Detection  of  Phases  by  Least  Squares  Filtering 

The  initial  P-raotion  is  probably  the  most 
easily  detectable  phase  on  a  seismogram.  The  signal  and 
the  microseismic  nois  in  general; posess  different  sta¬ 
tistical  properties  and  various  techniques  car.  be  used 
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to  enhance  the  signal-to-noise  ratio.  The  later  P-phases 
(such  as  pP,  PP,  PcP)  are  in  general  less  suitable  for 
purely  statistical  treatment ,  However,  these  later  events 
may  be  expressed  as  filtered  versions  of  the  initial  P- 
wave. 

Here  we  describe  a  scheme  of  detecting  these 
late  P-pirses  given  the  initial  P-wave  motion.  The 
technique  involves  finding  short  filters  which  match 
the  P-wave  to  later  phases  of  the  seismograms  by  mini¬ 
mizing  the  squared  error.  The  mean  square  error  determines 
the  quality  of  the  fit  and  can  be  used  as  the  detection 
criterion. 


The  simplest  formulation  for  the  noise  free 
case  can  be  given  as  follows: 


Let 

.  >  ^  i  ^ 1  ;  •  •  •  •  ^  a  ^ 

be  the  input  signal 

(i.  e.  initial  P-wave), 

be  the  filter  to  be 

determined  and, 

be  the  desired  output.  This  could  be  any  portion  of  the 
seismogram  where  we  wish  to  test  if  it  matches  the  initial 
P-wave.  Let  us  define  a  matrix 


r*.  ° . ° 

i  ** 
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Let  us  find  such  that  the  error 

£=  Xf-d 


has  minimum  power.  The  minimization  criteria  enable  us 
to  find  -j  .  The  mismatch  between  )(  and  Ji  can  be 
measured  in  terms  of  normalized  mean  square  error 


For  perfect  match  £  ~G  and  for  no  match  £  -  \  . 

In  practical  application  we  take  the  initial 
P-wave,  specify  the  length  of  j  ,  and  slide  it  along 
the  seismogram  at  each  point  computing  -£ •  and  £  . 

Where  there  is  a  matching  P-phase;  £  becomes  small. 

An  example  of  this  scheme  is  shown  as  applied 
to  LASA  recordings  (Figure  1A)  of  the  Long  Shot  event. 

In  Figure  15.  we  show  the  "performance  parameter" 

•  F  =--  f(Xrd  ) 

The  maxima  of  F  correspond  to  arrival  of  phases.  The 
crustal  reverberations  and  the  PcP  are  clearly  visible 
in  the  Figure  1A,  while  the  original  seismogram  does  not 
show  the  reverberations  as  clearly. 
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D. 


Excitation  of  Free 
bx  a  Point  Source 


—Oscillations  and  Surface  Wavpc 
— n  a  Vertically  Heterogeneous  Earth. 


Radiation  patterns  of  surface  waves  and  free 
oscillations  for  vertically  heterogeneous  elastic  media 
are  derived  for  arbitrary  sources  using  variational  equa¬ 
tions.  The  results  are  expressed  in  terms  of  normal  mode 
solutions  and  source  functions,  and  show  that  for  a 
spherical  earth  the  radiation  patterns  can  be  constructed 
from  normal  mode  solutions  alone.  Source  functions  for 
a  single  force,  a  single  couple,  and  double  ccuple  without 
torque,  all  in  arbitrary  directions,  are  derived. 
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IV.  SOURCE  MECHANISM  STUDIES 

A.  Scaling  Law  of  Seismic  Spectrum 

The  dependence  of  the  amplitude  spectrum  of 
seismic  waves  on  source  size  is  investigated  on  the  basis 
of  two  dislocation  models  of  an  earthquake  source.  One 
of  the  models  is  due  to  Haskell  and  is  called  the  OJ* 
model.  The  other  is  called  the  (jJ  model,  and  is  con¬ 
structed  by  fitting  an  exponentially  decaying  function 
to  the  autocorrelation  function  of  the  dislocation  vel¬ 
ocity.  The  number  of  source  parameters  is  reduced  to 
one  by  the  assumption  of  similarity.  We  found  that  the 
most  convenient  parameter  for  our  purpose  is  the  mag¬ 
nitude,  M  ,  defined  for  surface  waves  with  period  of 
s 

20  seconds.  Spectral  density  curves  are  determined  for 
given  Mg . 

Comparison  of  the  theoretical  curves  with  ob¬ 
servations  is  made  in  two  different  ways.  The  observed 
ratio  of  the  spectra  of  seismic  waves  with  the  same  propa¬ 
gation  path  but  from  earthquakes  of  different  sizes  are 
compared  with  the  corresponding  theoretical  ratio,  thereby 
eliminating  the  effect  of  propagation  on  the  spectrum. 

The  other  method  is  to  check  the  theory  with  the  empirical 
relation  between  different  magnitude  scales  defined  for 
different  waves  at  different  periods.  The  CO  model 
gives  a  satisfactory  agreement  with  such  observations  on 
the  assumption  of  similarity,  but  the  CO  J  model  does  not. 
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(See  Figure  16.) 

The  assumption  of  similarity  implies  a  constant 
stress  drop  independent  cf  source  size.  The  stress  drop 
in  earthquakes  may  h_  estimated  from  the  spectral  density 
of  surface  waves  at  long  periods,  if  the  fault  length 
and  offset  are  known  from  field  measurements .  A  pre¬ 
liminary  study  of  Love  waves  from  the  Parkfield  earth- 
equake  of  June  28,  1966,  shows  that  the  strtjs  drop  at 
the  source  of  this  earthquake  is  lower  than  the  normal 
value  (around  100  bars)  by  about  two  orders  of  magnitude. 
This  result  implies  that  a  single  parameter,  such  as  mag¬ 
nitude,  cannot  describe  an  earthquake  even  as  a  rough 
measure . 

Reference 
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B.  Earthquake  Prediction 

Recent  evidence  from  the  field  and  laboratory 
investigations  are  reviewed  to  make  a  critical  evaluation 
of  the  feasibility  of  earthquake  prediction.  Possible 
mechanisms  of  the  shallow  earthquakes  are  stated  and  the 
eifects  of  pore  pressure,  water  weaking,  lubrication,  and 
creep  fracture  are  discussed. 

There  are  some  field  and  laboratory  data  on 
the  premonitory  indications  of  an  earthquake.  Forewarning 
may  come  from  three  sources: 

1.  Observed  deformation  before  an  earthquake 
such  as  tilts  and  strain  in  epicentral  region  (These  have 
been  observed  and  documented  in  Japan). 

2.  The  general  increase  in  the  number  of  small 
events  prior  to  the  main  shock  (or  failure)  have  been  ob¬ 
served  in  the  laboratory. 

3.  Changes  in  physical  properties  of  rocks  near 
the  fault  may  take  place  as  they  are  strained.  In  the 
laboratory  a  drop  In  seismic  velocities  has  been  detected 
as  the  rock  stressed  to  the  point  of  failure.  Also,  a 
significant  decrease  in  the  electrical  resistivity  has 
been  measured  as  the  stress  is  increased  to  the  critical 
value . 


Suggestions  have  been  made  for  the  instrumentation 
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oC  seismic  belts  along  major  faults.  Type  of  instruments, 
and  field  clusters  and  general  field  investigations  have 
been  evaluated. 

Reference 

Press,  F.  and  W.  F.  Brace,  Earthquake  Prediction.  Science. 
152.  1575-1584,  June  17,  1966. 
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V.  SEISMIC  STATION 

A  seismic  station  to  be  equipped  with  long  period 
instruments  is  being  completed  at  Harvard,  Massachusetts. 
This  will  be  a  reference  station  for  long  period  record¬ 
ings,  and  testing  and  calibration  of  the  similar  field 
instruments.  The  seismic  vault  is  being  made  available 
by  the  Harvard  University,  At  the  present  some  repairs 
and  viodifications  of  the  vault  are  being  completeci. 

Initially  the  station  will  contain  a  two -cor.  - 
ponent  Benioff  tiltmeter  (Mercury  pendulum),  two -component 
quartz  rod  strain  meters,  and  a  long  period  vertical 
seismometer  with  displacement  transputer.  In  addition  to 
these  some  experimental  tilt  and  strain  devices  will  be 
placed  at  the  vault  for  testingc 

The  construction  of  the  tiltmeters  and  all 
analog  recording  systems  have  been  completed.  Currently 
a  number  of  digital  recording  devices  are  being  evalu¬ 
ated.  It  is  anticipated  that  analog  data  from  tiltmeters 
will  be  available  in  six  weeks,  and  those  from  strain- 
meters  in  less  than  three  months. 
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FIGURE  CAPTIONS 


Fig.  1 

Fig.  2 

Fig  3 

Fig.  4 

Fig.  5 

Fig.  6 
Fig.  7 

Fig.  8 

Fig.  9 

Fig.  10 

Fig.  11 


Locations  of  the  LASA  and  the  Long  Shot  nuclear 
explosion.  The  dashed  lines  outline  the  aximuth 
range  300° -3 20°  from  USA. 

Observed  dt/dA  for  all  events  in  the  azimuth  range 
300°-3?0c  detected  at  USA.  The  dashed  curve  is 
the  adopted  "mean"  curve. 

The  comparison  of  the  observed  dt/dA  ("good" 
events  from  Figure  7)  and  the  theoretical  values 
from  Model  77. 

Compressional  velocity  profile  of  the  Mantle  Model 
77,  with  the  Jeffreys -Bullen  model  for  compari¬ 
son.  Boxes  show  the  regions  of  the  "discontinui¬ 
ties"  . 

Observed  dt/dA  for  all  events  from  the  azimuth 
ranges  300-320°  and  140-160°  from  the  USA.  These 
deomonstrate  the  differences  between  the  two  dif¬ 
ferent  paths. 

Observed  composite  path  Love  wave  data  and  the 
extracted  "pure-path"  phase  velocities. 

Upper  mantle  shear  velocity  models  for  the  oceanic, 
continental  shield,  and  tectonic  regions.  Below 
500  km  profiles  are  uncertain  because  of  insuf¬ 
ficient  data. 

Temperature  differences  under  oceans  and  conti¬ 
nents  with  dV/dT  as  a  parameter.  A  constant  value 
(C)  can  be  added  to  or  subtracted  from  each  tem¬ 
perature  profile. 

Phase  and  group  velocity  dispersion  curves  for 
several  modes  oT  ARDC  standard  atmosphere  with 
free  surface  at  220  km. 

Phase  and  group  velocity  t  jpersion  curves  for 
fundamental  gravity  wave  in  a  constant-velocity, 
constant -density  ocean  of  5  km  thickness. 

Phase  velocity  dispersion  curves  for  modes  of  the 
air-ocean  system. 


Synthetic  barogram  (upper;  end  marigra^  (lower) 
for  San  Francisco.  Source  time  function  is  a 
single -cycle  sine  wave  of  40  minute  period.  Time 
is  local  civil  time,  August  27,  1883. 

Spectra  of  main  shock  (crosses)  and  main  plus 
aftershock  sequence  (circles)  frr  Chilean  earth¬ 
quake  of  1960.  For  this  case,  Q  assumed  at  200, 

Tp  *  10.0  minutes.  Peak  frequency  shift  is  .01 
minutes,  Q  shifts  to  250. 

Long  Shot  seismograms  of  A«  and  B  ring  of  LASA. 
Each  trace  is  a  composite  (sum)  of  25  seismometer 
outputs . 

"Performance  factor"  for  traces  shown  in  Fig.  14. 

Theoretical  relation  between  Mg  and  Mg  based  on 
the  W-square  model,  as  compared  with  that  observed 
by  Gutenberg  and  Richter. 
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